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Abstract. This paper is concerned with the effect of parameter fluctuations with 
a characteristic waiting time in coupled Rossler oscillators. We show that the 
averaged error in synchronization that is introduced due to a fluctuating parameter is 
proportional to the waiting time and the amplitude of the fluctuations. It is also shown 
that coupling strength beyond a threshold value does not have any significant effect 
on improving the quality of synchronization when the fluctuations posses considerable 
waiting time. 
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Synchronization of coupled chaotic systems has generated a lot of research activities 
over the last several years [USE]. Synchronized chaotic behavior has been studied 
extensively in physical, chemical and biological systems. One of the frequently employed 
methods is coupling two identical systems, this coupling may be unidirectional or 
bidirectional. Complete synchronization of identical chaotic systems is of considerable 
interest because of its applications in secure communication [U [6j [7J. By identical 
systems we mean a set of systems whose parameters are exactly equal. It is found that 
complete synchronization is not possible when there is a small but finite mismatch 
of the parameters of the systems. In coupled non autonomous systems effect of a 
phase mismatch or finite constant frequency detuning is to destroy the synchronization 
altogether [9|. 

Coupled Rossler oscillators have been extensively used for studies in synchroniza- 
tion. This is due to their simplicity and results obtained can usually be generalized to 
other chaotic systems. Kurths. et. al have shown that coupled Rossler systems which 
possess parameter mismatch attain a state of phase synchronization [TO] and on increas- 
ing the coupling strength, such systems may synchronize, but lagged in time [UJ. It is 
also found that many systems, both mathematical models and actual physical systems, 
show similar behavior [12J. There remains a question as to what happen if the param- 
eter mismatch is fluctuating. Such a study is relevant and important since in reality it 
is difficult, if not impossible, to construct identical systems except in numerical simu- 
lations. This can also be due to the fact that the parameters could be fluctuating in 
time with time scales of their own, either due to internal instabilities or due to some 
external perturbations. In this paper we discuss the effect of fluctuating parameters 
and the effect of timescales associated with such fluctuations on the synchronization of 
coupled Rossler oscillators. 

1. Parameter fluctuations 

Let us consider the parameters p\ and p 2 of a coupled system of oscillators which 
fluctuates randomly. The fluctuations are assumed to occur in time as follows 



where, £u and &t are two delta correlated zero mean random variables and po be some 
average value of the parameters. It is also assumed that the the fluctuations have a 
waiting time r . That is To is the time for which the parameter retains its value before 
it is changed once again. We define the fluctuation rate 0, as the number of times a 
parameter is perturbed in unit time which is also the inverse of the waiting time . In a 
case where the parameter fluctuates in this fashion, it can be seen the the parameter is 
correlated in time as 



Pit = Po + fit 
P2t =Po + £m 



(1) 
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Where p(..) denotes the normalized auto correlation. The amplitude of the parameter 
mismatch is denoted as Ap, given by 

Ap=<|fc(f)|>«, (2) 

where, 5p t = Pu—P2t and (...) t denotes time average. We did not consider the root mean 
square value because the system is sensitive to the magnitude of parameter mismatch, 
not to its square, at least in general. This model of fluctuation is, not the typical, but 
we believe that this is the closest approximation that can be applied an actual system. 
Also, it is interesting to note that our model of fluctuation is complimentary to the 
dichotomic noise [13] where amplitude of fluctuation is a constant, but the waiting time 
fluctuates about an average. 

There is one question that may arise at this point. The systems under consideration 
is continuous but the fluctuations are discrete. Why continuous fluctuations, for 
example, an Ornstein Uhlenbeck kind, which posses time scales and are continuous 
were not considered. Our answer is, such fluctuations arise as a result of a stochastic 
evolution. As far as complete synchronization is considered, we do not expect such 
untethered stochastic evolution of parameters in the parameter space. Instead, a 
parameter, if deviated from its desired value is expected to fluctuate around an average 
as complete synchronization is usually not found in nature but found in fabricated 
systems, well monitored and well designed. There, we expect if a parameter mismatch 
occurs, it will have a tendency or it will be forced to cross zero very often. In 
such a situation the prime consideration shall be the time scales associated with the 
fluctuations. 



2. Effect of Fluctuations on Dynamics 

In a case where the parameter fluctuating as described in sec. Q] the effect on the 
dynamics of coupled chaotic system can be understood in terms of the dynamical 
equations. Let the evolution of the coupled systems in phase space be given by, 

X 1 = f 1 (p 1 ,X 1 ) + Cf(X 2 -X 1 ) (3) 
X 2 = f 1 (p 2 ,X 2 ) + Cf(X 1 -X 2 ). 

Where X represents the phase space variables, p the parameter whose fluctuation is 
considered, and C, the coupling strength. With equation (jBJ) we can write an equation 
for the rate of separation X\ — X 2 of the trajectories as, 

d{Xl ~ t X2) =X,-X 2 = M(p 1 ,p 2 , X 1 , X 2 ), (4) 

M(pi,p 2 , Xi, X 2 ) is a function of the dynamical variables, the parameters of the coupled 
systems and Ap the parameter mismatch. This can be expanded int terms of Ap and 
the effect of fluctuations can be separated out. 

M(p 1 ,p 2 ,X 1 ,X 2 ) = M s (p ,X u X 2 ) (5) 
+ E(X 1 ,X 2 ,Ap 1 ,Ap 2 ). 
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Here M s (p , Xi, X 2 ) represents the quantity which offers a stable synchronization 
manifold and E(Xi, X 2 , Api, Ap 2 ) represents the effect of the parameter mismatch. 

Let us now see the form of E(Xi, X 2 , Api, Ap 2 ) in the in case of a system of 
bidirectionally coupled Rossler oscillators. 

xi,2 = - 2/1,2 - zi,2 + c(x 2 ,i - x 1>2 ) (6) 

Vl,2 = %1,2 +Pl,2Vl,2 

z li2 = 0.2 + z li2 (x 1;2 - 10) 

Here it is assumed that in the absence of fluctuations and for an appropriate value of 
c, the systems get synchronized. Also in the presence of fluctuations an approximate 
synchrony is maintained due to the negative conditional Lyapunov exponents and the 
zero mean nature of the fluctuations. With equation ([6]), we can write the rate of 
separation of trajectories as, 
d(x 1 - x 2 ) 



dt 

d(yi - 2/2) 
dt 

d{z x - z 2 ) 
dt 



- 2/1 + 2/2 - zi + z 2 + 2C(x 2 - xi) (7) 
Xi - x 2 + p x yi - p 2 y 2 
X\Z\ - x 2 z 2 - 10(2:1 - z 2 ). 



Here it can be seen that the fluctuations affect the dynamics through y, as the parameter 
p occur only in the equation for y in the coupled set of equations. Thus by Taylors 
expansion around po, E(X\, X 2 , Ap\, Ap 2 ) can be written as, 

E{X X , X 2 , A Pl ,Ap 2 ) = y^ x - y 2 & (8) 

The evolution of the system in time can be split into several sub intervals of duration 
t. Thus in a approximately synchronized state the error occurred due to the fluctuations 
in the k th interval can be written as, 

rt+T 

e{t,t + r)= jf {yxix-y^dt (9) 

If the parameter fluctuation are fast enough yi and y 2 can be considered to be a constant 
between two consecutive fluctuations and is denoted by yi and 2/2- Also £1 and £1 are 
constants for a given waiting time r by definition, and their value in the k th interval is 
and £ 2 jk- Thus equation ([9]) can be written as, 

rt+T 

e k (t,t + T )= j (y!k£ik ~ y2k&k)dt (10) 
= To(2/Tfe£ifc - 2/2fc6fc) 
From this equation, (e^) can be calculated 

(4) = r 2 (m 2 e lk + m 2 e2 k ) (ii) 
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In the above relation (yTtf^fk + V^Qk) can be considered as a constant as y belong 
to an attractor which is dense in periodic orbits and the averages of £ over the waiting 
time follows some well defined statistical distribution. (e\)t can be written as 

(el) = Tq x constant (12) 

as t — > oo. Here note that this expression do not contain C which is the coupling 
strength. This suggest that increasing the coupling strength do not have any significant 
effect on reducing the fluctuations. Also note that the error introduced by fluctuations 
is proportional to the amplitude of fluctuations. 



3. Effect of Fluctuations on the Quality of Synchronization 

In the last section we saw that the the fluctuations affect the dynamics in a multiplicative 
manner and its effect on the seperation of the trajectories is proportional to the waiting 
time. In this section we consider how the fluctuations affect the over all dynamics of the 
coupled system. Consider a system where the dynamics is represented by the dynamical 
equations z>j = f i (i>i...v n ). The effect of a perturbation of a variable Vf. on the variable Vj 
can be written as 5uj = ^j^f^- Sv^. Thus it can be seen that in a situation where the 
phase space variables can be considered to be a constant, the effect of the perturbation 
applied to one variable have proportional effect on other variables also. 

Effect of such perturbations on the quality of synchronization can be quantified 
using any measure of synchronization which is based on the divergence of the trajectories 
of the coupled system. A well known measure of the quality of synchronization which 
is the similarity function 5(0) defined as, 

S 2 (Q) = (Mt)~X 2 (t)] 2 ) 
[<S?(t)><*2(t)>]* " 

In an ideally synchronized state this directly corresponds to e^. If an approximate 
synchrony is maintained during the evolution of the system we can assume that the 
majority of contribution to the error in synchronization is of the form e^. Thus the 
analysis that we have presented is valid if 5(0) can be fitted with a function which 
is linear in r , or proportional to the inverse of cj), the fluctuation rate in an actual 
numerical experiment. 



3.1. Numerical simulations 

In figure [1] it can be seen that the coupling strength greater than the threshold value 
does not play any significant role in determining the quality of synchronization. But 
the quality increase as the fluctuation rate is increased. In figure [2] curve fitting is done 
for the coupling strength c = 0.9, it can be seen that 5(0) varies with To as 

S(0) = a + br . (14) 

where a = 0.005 and b = 0.12 In figure [3] it is shown that for two fluctuation rates, 5(0) 
grows linearly with the amplitude of fluctuations, the growth rate of the error is higher 
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Figure 1. The synchronization error decreases with the increase in the fluctuation 
rate. It can be seen that high coupling could not stabilize synchronization with lower 
fluctuation rates. Here Ap = 0.05. 



for a larger waiting time. Numerically, the results were similar for the averages of £ over 
the waiting time following uniform or Gaussian distribution. Thus for a given amplitude 
of fluctuations it is the time scales associated with the fluctuations that determine the 
quality of synchronization. Interestingly there are similar results in biological systems 
[15] with coloured fluctuations that higher correlation times (time scales) makes the 
coupled systems less synchronizable. 

With a low fluctuation rate, the parameter fluctuations can considerably affect 
synchronization because the phase space evolution time is comparable to the interval 
where a fixed parameter mismatch persists. Thus, the system always get time to respond 
to the parameter mismatch before it being canceled out. The error-timescale relations 
in this regime can be different, which requires further work. 

4. Conclusion 

In this paper, we show that the synchronization error introduced due to a deviation of the 
parameter from its desired value, is proportional to the waiting time, the instantaneous 
value of the variables and the amplitude of perturbation. Asymptotically this leads to 
a relation between the similarity function and the fluctuation rates which is reciprocal 
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Figure 2. Relation between <j> and 5(0) is found to be of the form 5(0) = a + b * t , 
a=0.005 and b= 0.12. 



in nature. Also the coupling strength which plays an important role in determining 
the nature of synchronization in systems with constant parameter mismatch, but do not 
have any significant role when the mutual parameter mismatch is fluctuating. It is hoped 
that this investigation will spur further research in this field, from a more fundamental 
point of view as well as for practical implications where high quality synchronization is 
required. 
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Figure 3. For two fluctuation rates, </> = 50 (red) and <j) = 500(black): the 
synchronization error grows as the amplitude of fluctuation is increased. Note the less 
significant role of coupling strengths compared to that of fluctuation rates. 
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